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Abstract. In this paper we present the WeCAPP variable star catalogue towards the bulge of M31. The observations 
of the WeCAPP microlensing survey (optical R and / bands) during three years (2000-2003) result in a database 
with unprecedented time coverage for an extragalactic variable star study. We detect 23781 variable sources in 
a 16.1' X 16.6' field centered on the nucleus of M31. The catalogue of variable stars contains the positions, the 
periods and the variations in the R and / bands. We classify the variables according to their position in the _R-band 
period-amplitude plane. Three groups can be distinguished; while the first two groups can be mainly associated 
with Cepheid-like variables (population I Cepheids in group I, type II Cepheids and RV Tauri stars in group II), 
the third one consists of Long Period Variables (LPVs). We detect 37 RV Tauri stars and 11 RV Tauri candidates 
which is one of the largest collections of this class of stars to date. The classification scheme is supported by 
Fourier decomposition o f the light curves. O ur data shows a correlation of the low-order Fourier coefficients "l>2i 
with <l>3i (as defined bv lSimon fc Leel (|l98l')) for classical Cepheids, as well for t ype II Cephei ds a nd RV Tauri 
stars. Correlating our sample of variable stars with the X-ray based catalogues of lKaare'3 i2002l) and lKong et alJ 
||20o2) results in 23 and 31 coincidences, 8 of which are M31 globular clusters. The number density of detected 
variables is clearly not symmetric, which has to be included in the calculations of the expected microlensing event 
rate towards M31. This asymmetry is due to the enhanced extinction in the spiral arms superimposed on the bulge 
of M31 which reduces the number of sources to about 60%, if compared to areas of equivalent bulge brightness 
(without enhanced extinction being present). 
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1. Introduction 

In the last decade the ongoing microlensing surveys 
greatly extended our knowledge of variable stars. As 
the observations of these experiments usually cover a 
long time span with a good time sampling the resulting 
datasets are suited perfectly for the study of many differ- 
ent types of variable sources. Lots of the progress in this 
field therefore result ed from the work o f collaborations 
like MACHO (e.g., lAlcock et~al] lll995[ [1998)), EROS 
(e.g., iBeauHeu et al." ("1995"); 'Derue et al.' ^2 



(e.g., ICieslinski et al. 



Derue et al.' (2001 ,)), OGLE 
(2001) : Wrav. Ever fc Paczvnskil 
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t()Q4\ ) or MOA (e.g., iNoda et alJ ll2002l l2004l) 1. 

Numerous publications enlarged the list of known vari- 
able stars like Cepheids, RR Lyrae or galactic Long Period 
Variables (LPVs), but also helped to understand the phys- 
ical processes dominating these stars. 

M31 was surveyed for variable sources since the 1920s 
starting with the pioneering work of Edwin Hubble. With 
plates taken at the newly available Mount Wilson tele- 
scope Hubble succeeded to resolve Cepheid variables in the 
outer parts of M31. Using the alre ady established period- 
luminosity relation, iHubbld l)l922j) was able to determine 
the distance to M31 to 300 kpc (the difference to the actual 
value of 780 kpc (Sta nek fc Garnav ich. 1998) is mainly 
due to an erroneous calibration of the zero-point of the PL- 
relation and a missing reddening correction) and in this 
way to reveal the extragalactic nature of the Andromeda 
Ne bula'. 

iBaade fc Swopd lll96.1 ll96Rl) continued the work on 
variable stars in M31 and detected over 400 variables, 
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among them Cepheids and novae. Looking at the rela- 
tions between period, luminosity, amplitude and frequency 
they found that the Cepheids resemble the ones in the 
Milky Way, but seem to be different from those in the 
Small Magellanic Cloud. Another comp rehensive study 
was p erformed by the DIRECT project ijKaluznv et al.1 
Il998l) who detected and examined Cepheids and Detached 
Eclipsing Binaries in five fields towards the M31 disk in 
order to reduce the uncertainties in the distance deter- 
mination of M31. Unfortunately, the data of a 6th field 
towards the bulge of M31, which in part overlaps with our 
surveyed area, remained unpublished. 

In the last years Andromeda wa s target of 
sever al microlensing su rveys (AGAPE llAnsari et all 
I1997I) POINT - AGAP E llAuriere et all l200lh. W eCAPP 
llRiffeser et aD. 120011). MEGA llde Jong et all |2004), 
SLOTT/AGAPE l|Calchi Novati et al.ll2002t) ) looking for 
compact dark matter in the halo of M31. As by-product 
of these surveys, catalogues of variable sources start to 
appear. The work presented in this paper overlaps with 
the work resulting froi n the AGAPE Ansari et aill2004|) 
and POINT- AGAPE l)An et al.Ll2004 data sets. 

Due to our daily stacking scheme of the observations, 
we are not sensitive to variations and periods smaller than 
1.3 days and focus therefore on variations of longer pe- 
riod variables, namely Cepheids, RV Tauri stars, and Long 
Period Variables (LPVs), i.e. Miras and semi-regular vari- 
ables. The catalogue is completed by eclipsing binaries, 
and variables showing eruptive or irregular variations. S 
Scuti and RR Lyrae stars show variation amplitudes be- 
low our detection limit and will therefore be missed in this 
study. 

The paper is organized as follows. In Sec. |21 we give 
an overview of the survey, the observations, and data re- 
duction. Sec. 121 deals with the source detection and the 
derivation of possible periods of the variables. In Sees. 01 
and El we show the different groups of variable sources 
detected in the survey. The catalogue of variable stars is 
presented in Sec. (HI and is correlated with X-ray selected 
catalogues in Sec. (71 Whereas Sec. |H1 summarizes the pa- 
per, the appendix deals with the accuracy of the derived 
periods. 



2. The dataset 

2.1. Description of the survey 

The data presented here result from three years observa- 
tions of the central part of M31 by the WeCAPP project. 
We obtained data from the 0.8 m telescope at Wendelstein 
Observatory (Germany), and from the 1.23 m telescope at 
Calar Alto Observatory (Spain). At Wendelstein with its 
field of view (FOV) of 8.3' x 8.3' we selected a field (Fl in 
the following) along the minor axis of M31 which contains 
the area with the highest expected rate for pixellensing 
events (see also Fig.^for the location of Fl). Observations 
of this field were accompanied by images of F3, the oppo- 
site field along the NW minor axis, taken with a sparser 



Fig. 1. V-,R-, and /-band composite image of the ob- 
served fields Fl to F4, taken at Calar Alto Observatory 
during the campaign 2000/2001. The black lines mark 
the positions of fields Fl to F4. The identified variable 
sources lie within RA(2000): [00h43m25.0s, 00h41m59.9s] 
and DEC(2000): [41d08'00.1", 41d24'18.0"]. 



time sampling. The Calar Alto field covered 17.2' x 17.2' 
centered on the nucleus of M31. Two quadrants of the field 
coincided with the maximal lensing field Fl and the op- 
posite field F3. Due to the simultaneous observations we 
reached a very good time sampling during the observabil- 
ity of M31. Since summer 2002 we are mosaicing fields Fl 
to F4 with the Wendelstein telescope solely. A composite 
image {V-, R-, and /-band) of fields Fl to F4 taken at 
Calar Alto Observatory during the campaign 2000/2001 
is shown in Fig. ^ The epochs with data taken for the 
four fields are shown in Fig. [3 
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Fig. 2. Distribution of the observations for the four 
fields. During the first two campaigns we reached a very 
dense time sampling of the observations. Note that the 
third and fourth campaign were restricted to Wendelstein 
Observatory only. 
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2.2. Instruments and filters 



An overview of the cameras and CCDs used in WeCAPP is 
shown in Tabled The cameras are a IK x IK TEK cam- 
era with a pixel size of 0.49" at Wendelstein and a 2K x 2K 
SITe camera with a pixel size of 0.5" at Calar Alto. The 
observations were carried out in R and / filters close to the 
Kron-Cousins system. At Wendelstein we used the i? (A ~ 
650 nm, AA ~ 150 nm) and / (A ~ 850 nm, AA ~ 150 nm) 
wavebands. The Calar Alto observations were carried out 
with the very similar filters, R{X ~ 640 nm, AA ~ 160 nm) 
and / (A ~ 850 nm, AA ~ 150 nm). Despite of the com- 
bination of different telescopes, CCDs, and slightly dif- 
ferent filter systems we observed no systematic effects in 
the light curves depending on these parameters (see also 
Fig. [2H1 where we show the data obtained at Calar Alto 
and at Wendelstein separately.) 
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2.3. Data reduction 

The WeCAPP reduction pipeli ne mupipe will be pre- 
sented in a f uture pubhcation i Riffeser. Fliri. fc GosslL 



12005. in prep.h . mupipe combines all reduction steps from 
de-biasing of the images until PSF (point-spread-function) 
photometry and the final measurements of the light-curve 
in one software package, including full error propagation 
from the first reduction step to the last ijCossl &: Riffeseii 
I2OO2I) : 



j dttfereoce frame iinai error fram& 



PSF-f Stfing -phcKofneay 



1. standard CCD reduction including de-biasing, flat- 
fielding and filtering of cosmic ray events 

2. position alignment on a reference grid using a flux con- 
serving interpolation routine 

3. photometric alignment of the images 

4. restoration of pixels damaged by CCD defects (cold 
and hot pixels, traps, bad columns) or pixels hit by 
cosmic ray events 

5. stacking of the frames of one epoch (i.e. night) using 
a weighting scheme to maximize the S/N for point 
sources 

6. matching of the PSF of a high S/N reference frame 
to the PSF of each stacked frame using our implemen- 
tation of the image sub trac tion method developed by 
lAlard fc Liiplm] ^ll998^ and lAlardI ^l2000^ 

7. generation of difference images by subtracting the con- 
volved reference frame from each stacked frame 

8. PSF photometry of each pixel in the difference frames 
using a PSF extracted from the convolved high S/N 
reference frame 

9. extraction of the light curves 

A floating chart of the reduction pipeline is shown in 
Fig. O In the last step mupipe returns roughly 4 xlO^ 
pixel light curves together with appropriate errors, each 
of the light curves representing the time variability of 
the flux present inside the PSF centered on the particu- 
lar pixel. The extraction of intrinsic variable sources from 
these pixel light curves is presented in the next section. 



Fig. 3. Floating chart of the WeCAPP reduction pipeline 
mupipe. The reduction package includes full error propa- 
gation for each pixel through all reduction steps. In this 
way, all data points are properly taken into account in 
the search for variable sources. The reduction pipeline re- 
turns light curves for each pixel in the frame, representing 
the temporal change of the flux present inside the PSF 
centered on the particular pixel. 



2.4. Astrometry 

The astrometric solution was created using the IRAF^ 
tasks ccmap and cctran for 92 stars whose positi ons were 
taken from the Local Group Survey (LGS, Massev et alJ 
(2001)). The rms of the solution is 0.10" in declination and 
0.12" in right ascension. The errors of the astrometric so- 
lution for the 92 calibration stars are shown in Fig. 01 The 
coordin ates derived he r e agre e per fectly with the one's de- 
rived in lRiffeser efaH l^2001^ and lRiffeser et alJ l|2003() . 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 



4 



Fliri et al.: WeCAPP variable stars 



Site 


Campaign 


CCD 


Size 


[arcsec / px] 


Field [arcmin^] 


We 


2000-2003 


TEK #1 


IK X IK 


0.49 


8.3 X 8.3 


CA 


2000-2002 


SITe2b #17 


2K X 2K 


0.50 


17.2 X 17.2 



Table 1. Properties of the CCD cameras used during WeCAPP at Wendelstein (We) and Calar Alto (CA) 
Observatories, respectively. Both CCDs have a pixel size of 24 /im. 



ARA[arcsec] 



Fig. 4. Errors of the astrometric solution for the calibra- 
tion stars taken from the LGS survey. The rms is 0.10" in 
declination and 0.12" in right ascension. 



3. Detection of variable sources 

3.1. Selection of the sources 

For the selection of the intrinsic variable sources we use the 
i?-band data from the Calar Alto campaigns 2000/2001 
and 2001/2002 to create a xt variation frame. To get rid 
of systematic effects induced by the different seeing con- 
ditions of each frame we use the following approach. We 
select a stacked image with roughly the median seeing of 
the 2000/2001 campaign (about 1.5") as master frame and 
match the PSF of each stacked frame with a smaller PSF 
to the one of the master frame. In this way we end with 
difference images for 113 epochs (i.e. 48 % of the epochs 
with i?-band data obtained at Calar Alto) which are used 
for the detection of the variable sources. For each of the 
light curves extracted from the PSF1.5 difference images 
we calculate the reduced xt deviation from a constant 
(i.e. zero) baseline fit. The errors entering this calcula- 
tions are the propagated errors as returned from mupipe. 
The results are written into a x^-frame of the field. The 
mode of the x^-frame is 1.02 which shows the accuracy 
of the propagated errors (see Fig. [SI). Each value > 1.15 
is connected to a non-constant, i.e. variable source at the 
99.99 % confidence level. 

To avoid contamination of the M31 sample with fore- 
ground objects we set all pixels in a radius of 5 pixels 



around the positions of bright foreground stars to zero 
before detecting the variable sources in the x^-frame. We 
determine the positi ons of the variable source s using the 
SExtractor software ijBertin fc Arnout3 . Il996|) for source 
detection applied on the x^-frame. In this way we de- 
tect 25571 variable sources in our observed field. FigureEl 
shows an extract of the x^-frame with the positions of the 
detected sources marked. 



3.2. Period determination 



We use an a lgorithm developed by iLombl l)l976|) and 
IScargl to determine the significance and value of 

a possible period of the variables sources. The method of 
Lomb is well suited for this problem as it is able to deal 
with unevenly sampled data. 

The algorithm extracts the power in the first sine and 
cosine terms for a set of equidistantly spaced frequencies 
and also yields the significance of the detected peaks in 
the power spectrum. The Lomb normalized periodogram 
for N measured data points hj — h{tj), j = 1, N, taken 
at epochs tj is defined by 



Y,Ahj - h) cosaj{tj - t) 



1 2 



2cr2 



cos^ uj {tj 



h) sin Lo {tj 



2a2 



shi UJ {tj - r) 



(1) 




pixel value 

Fig. 5. Histogram of the pixel values measured in the x^- 
frame. The maximum of the distribution of 1.02 refiects 
the accuracy of the propagated errors. 
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Fig. 6. Extract (1' x 1') of the x^frame of the over 25000 
variable sources detected by the WeCAPP microlensing 
survey of M31. The circles give the positions of the sources 
detected by SExtractor. Most of the sources are first 
identifications, the majority of them being Long Period 
Variables (LPVs), i.e. Miras and semi-regular variables 
(see Sec. El). 



where uo = 2'kv = ^ is the angular frequency for the 
period P, and the mean h, the variance cr^, and the con- 
stant r are defined as follows 



1 ^ 
h = — > hi , 



N 



1 ^ 



tan(2wr) 



sin 2ijjtj 
^ cos 2ujtj 



(2) 



(3) 



(4) 



The last definition ensures the power spectrum to be in- 
dependent on a shift of ea ch by A^, as t transforms 
in that case into (r + At) . iLombl l|l976fl has shown that 
the evaluation of the periodogram according to Eq. ^ is 
identical to a linear least-squares fit of the first harmonics 



F{t) — A sin uit + B cos uit 



(5) 



to the data points. In Fig.[7|we show two examples for 
P/v(w), the first one derived from a Cepheid hght curve, 
the second one derived from the light curve of a Long 
Periodic Variable. 

The significance level P{> z) of a peak with amplitude 
z in the power spectrum is given by 



P{> z) = 1 - (1 



(6) 



with (1 — e ^) being the probability that none of the 
M tested frequencies shows an amplitude greater than z 



3 
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Fig. 7. Two examples of power spectra Pn{(^) of 
WeCAPP variables. In the upper panel we show the pe- 
riodogram for a Cepheid (see also Fig. I24|l derived from 
the i?-band data, in the lower panel we show the peri- 
odogram for a Long Periodic Variable (see also Fig. I28|l 
derived from the /-band data. 



in case of pure Gaussian noise. iHorne fc Baliima^ ljl98fj) 
showed with extensive Monte Carlo simulations, that the 
value of M, i.e. the number of independent frequencies 
does not differ much from the number of data points A'^, 
if the data points are not closely clumped. 

We use the implementatio n of the algorithm taken 
from the Numerical Recipes llPress et al. . 1988 ) . Small 
modifications in the code allow us to search for the dif- 
ferent maxima in the power spectrum. The positions and 
values of the peaks are returned and can be used for fur- 
ther study of the light curves. In this implementation the 
significance P{> z) is approximated for small values of 
P{> z) by 



P{> z) w Me" 



(7) 



where M is a product of the number of data points N 
and a user supplied value, which determines the high end 
cut-off of the tested frequencies relative to the Nyquist 
frequency. 

4. The variable sources 

4.1. Creating the catalogue 

We detect 25571 variable point sources in a 16.1' x 16.6' 
field centered on the nucleus of M31. As this sample still 
can contain spurious detections or sources with only a 
few data points measured in the light curve we apply a 
couple of cuts which mark different levels of accuracy of 
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Fig. 8. Histogram of the number of epochs of the level 
A sources in the R— band (black curve) and /-band (red 
curve). For sources located in field Fl we usually got data 
points in the i?-band on more than 400 epochs. 

the derived periods and the required need for a visual 
inspection. 

As first cut we demand that the light curve comprises 
at least 40 data points in both filters. All sources which do 
not pass this criterion are removed from the sample. This 
reduces the number of variables under inspection to 25316 
and defines the level A sample. In Fig. |H1 we show the 
histogram of the number of epochs for these sources in the 
i?-band (black curve) and /-band (red curve). The lowest 
number peaks represent the time sampling present in fields 
F2 and F4, followed by the peak corresponding to field F3. 
As F2 and F4 were observed with comparable frequency 
both fields contribute to the lowest number peak. The last 
peak with usually more than 400 epochs in the /?-band is 
related to field Fl. 

As the Lomb algorithm assumes that all data points 
have the same error (see Eqs.^andl^l, data points with 
big error bars or outliers can spoil the period finding pro- 
cess and yield spurious periods. To avoid this we first elim- 
inate in each light curve ten data points with the biggest 
errors and in a second step five data points with the high- 
est and lowest values, respectively. In doing so, we ensure 
that we have at least 20 data points in both filters to look 
for a periodic signal. 

We now check the /?-band and /-band light curves 
separately for periodicities. As the color of variable stars 
usually changes during a cycle combing the two data sets 
would not result in accurate determinations of the period 
of the sources. Finally we obtain 25316 variables with de- 
termined periods (of any significance level) in both bands. 

We regard those periods as real which are the same 
in both bands inside tight boundaries reflecting the error 
of the period determination. For periods Pr < 215 days 
we take the theoretical error (see Eq. IA.2(I resulting from 
vpwHM as limit, from 215 days onwards we choose a more 
conservative constant limit of 30 days, which is increased 
to 60 days for Pr > 400 days. This criterion defines the 
level B sample (20311 objects). 

As final cut we select all variables from level B which 
show a significance level of the period determination of 
P{> z) < 10^^" in R or / and which have a determined 
/?-band period Pr < 450 days. This final cut reduces the 
number of sources to 19551 and defines the level C sam- 



ple with well determined periods. The remaining sources 
in level B were inspected visually and assigned to the 
group of regular, irregular, long- variation (period could 
not be determined because of an incomplete cycle) or mis- 
cellaneous variables in the final catalogue, or rejected 
from the sample as spurious detections. Finally, all sources 
of the level C sample are added automatically to the final 
catalogue. 

For the variables which show deviant periods in R and 
/ we proceed as follows. If the significance level of one 
of the deviant periods is better than 10~^^ and at the 
same time better than the significance in the other band 
by a factor of 10^° we choose the period in the first band 
and add the variable to the final catalogue. If this is not 
the case we inspect the light curve visually and decide 
if one of the periods is the real one, or if the variation 
is of irregular nature. In this step we also reject the last 
spurious detections from the sample, the true variables are 
added to the final catalogue. 

4.1.1. Search for eclipsing binaries 

By visual inspection of the folded light curves of the Lomb 
sample we detect 28 eclipsing binary (EB) candidates, 
amongst them one (semi-)detached system. However, the 
Lomb algorithm can fail in detecting periods for potential 
EB systems in M31. As EB light curves, especially the ones 
of detached systems, show strong power in the higher har- 
monics, they are not well recovered by Fourier techniques 
using the first harmonic only. Since some EB candidates 
therefore could be missed by using the Lomb algorithm 
alone, we add another step to the search for EBs. We 
investigate all light curves which were removed in the pre- 
vious steps once more, however use this ti me the dedi- 
cated transit finding algorithm 'boxfitting' ( Kqvacs et al.l 
I2OO2I) for the period determination. As lTinglevI l|2003() has 
shown the boxfitting algorithm is a powerful tool for the 
detection of eclipsing systems. We implemented the origi- 
nal Fortran code^ into our detection pipeline and run the 
algorithm over the mostly noisy light curves which were 
removed from the sample as spurious detections in the 
previous steps. In this way we assure to have proper pe- 
riods for all potential EB candidates present among our 
variable sources. 

Selecting all light cu rves with a reasonab le signal de- 
tection efficiency (SDE, iKovacs et aP ()2002j) ') SDE > 6, 
periods P > 1.30 days, and rejecting periods between 1.95 
days and 2.05 days (to avoid to pick up aliasing periods) 
results in an additional sample of 155 light curves. Visual 
inspection of the folded light curves finally yields 3 ad- 
ditional candidates for (semi-)detached eclipsing binary 
systems. 

In Fig. El we show the distribution of periods of the 31 
EB candidates in the upper panel. The 4 (semi-)detached 
systems populate the low period area up to a period of 7 
days whereas the candidate contact systems generally have 

^ http:/ /www.konkoly.hu/stafT/kovacs.html 
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Fig. 9. Top panel: Histogram of the periods of the 31 
echpsing binary candidates. The 4 (semi-)detached sys- 
tems populate the low period area up to a period of 7 
days whereas the candidate contact systems generally have 
longer periods. Bottom panel: binned _R-band light curve 
of a semi-detached system in the phase representation. 
The system has a i?-band period of 7.08 days. 



longer periods. In the bottom panel we show the binned 
i?-band light curve of a semi-detached system in the phase 
representation. The system was detected using the Lomb 
algorithm and has a i?-band period of 7.08 days. 

4.2. Number counts and asymmetry 

In Fig. ^1 we show the positions of the variables from 
the final catalogue, which suggests a connection of the en- 
hanced extinction in M31's spiral arms and the depletion 
of sources in certain regions in the northern part of the 
bulge (field F3). This depletion is also evident from num- 
ber counts of the sources in the northern and southern 
hemisphere of M31 (see Fig lll(l . 

We are calculating the number density of variable 
sources by summing up all variables in 100x100 pixel 
[50" X 50"] bins. In the resulting density map (Fig. [T^ 
the spiral arms are clearly visible. A comparison of the 
number densities at equivalent positions in the M31 bulge, 
one in the dust lane of the spiral arms, the other one in 
the opposite hemisphere, shows a reduction to about 60% 
relative to the part with no strong extinction. 

To examine this subject further we compare the po- 
sitions map (see Fig. IIOII with an extinction map of our 
field. In Fig. ^] we plot the positions of the sources on top 
of the i?-band extinction map, which we derived from the 
V- and i?-band frames taken by the LGS survey in the 
following way: 
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Fig. 10. Positions of the 23781 variable sources of the 
final catalogue. Blank regions are connected to to satu- 
rated parts of the frame or to eliminated foreground stars. 
Incompleteness induced by crowding in the center and the 
influence of the enhanced extinction connected with the 
spiral arms superimposed over the M31 bulge are clearly 
visible (see also Fig. 112(1 . 
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Fig. 11. Number counts of sources of the catalogue as 
function of the distance to the center of M31. Black: all 
sources. Blue: sources in the northern hemisphere of M31. 
Red: sources in the southern hemisphere of M31. Clearly 
visible is the asymmetry of the detected sources due to 
the enhanced disk-extinction in the northern part. The low 
number of sources near the center is due to incompleteness 
induced by crowding and saturated parts in the frames. 



We start from the relation between the color excess 
E{ V — R) and the extinction Aj^ in the i?-band (a « 3, 
see lBinnev .sT Mcr rifield (,1988,\ ) 



Ar = a EiV - R) 



(8) 



With the non-reddened magnitudes Mii,o,Mv,o and red- 
dened counterparts Ma^r, Mv,r this writes as 



M 



R.r 



Mr^O = a [iMv,r - Mr,.) - (Afy,o " Mfl.o)] • (9) 



8 



Fliri et al.: WeCAPP variable stars 




Fig. 12. Number density of variable sources by summing 
up all variables in 100 x 100 pixel [50" x 50"] bins. To avoid 
effects induced by the border of the frame we only use 
pixels in the region [200:2000,200:2000]. The contour lev- 
els are 1, 2, 3, 4 and 5 xlO~^ [arcsec"^], respectively. The 
spiral arm is clearly visible. Due to extinction in the dust 
lanes the number density reduces to about 60% when com- 
pared to equivalent regions without strong dust features 
being present. The regions with underdensities near the 
center are due to saturated parts in the frame with no 
variable sources detected. 



The transformation to fluxes Fij [i — R, V; j = 0, r) 
yields 



Frh = F] 



R,r 



FR.r\°' ( Fyx 



'V,r 



(10) 



If we now suppose that the intrinsic stellar popula- 
tion gradients are negligible over the field we can set 
('^^) ^ constant, which is not exactly true but a valid 
approximation for our purposes. 

Using this assumption we finally obtain a relation for 
Fro 



Fr,o — 



F. 



R,r 



pa 



Fr^ 



(11) 



which can be used to calculate the extinction extR in the 
i?-band 



extjj. = —2.5 log 



/ Fr^t 



\Frx 



(12) 



For the intrinsic color {V — R)q required in the calcula- 
tion we use the theoretical value {V — R)o = 0.63, assum- 
ing the bulge to be a 12 Gyr old SSP of 2Zq metalicity (C. 
Maraston, priv. comm.). Figure lTSI shows the resulting ex- 
tinction map: indeed under-dense variable source regions 
coincide with high extinction regions. 




Fig. 13. i?-band extinction map of the M31 field. White 
areas correspond to high extinction values. The extinc- 
tion in the JJ-band was derived using t he V- and it!-ban d 
images taken from the LGS survey (M assev et alil20oH) . 
Overplotted are the positions of the catalogue variables 
(blue dots). In the central regions the extinction map 
could not be calculated due to saturation in the original 
LGS frames. 



Finally we compare the number densities of sources 
present in our ca talogue with theoretical predictions. 
According to lli.enzini. (.1998.) the number of LPVs per inte- 
grated 10^ bolometric luminosities is equal to 0.5. lR.enz'in3 
1)19981) uses a slightly older (15 Gyrs) and less metal-rich 
(1 Zq) bulge if compared to the model used for the cal- 
culation of the e xtinct i on fra me. With the small bulge 
decomposition of iKentI lIlQS Q*) and under the assumption 
that the LPVs are present in the bulge only, we derive the 
expected number densities of LPVs and show it in Fig. 1141 
The agreement with the detected number densities is good 
in the outer parts of the field. Towards the center we suf- 
fer from incompleteness due to enhanced noise on the one 
hand and crowding of the sources on the other hand. 

4.3. Classification scheme 

As practically all of the detected sources are unresolved in 
the original frames we are neither able to put the sources 
in the color-magnitude plane and to construct a color- 
magnitude diagram, nor to derive period-luminosity rela- 
tions for the different classes. The available parameters 
for establishing a classification scheme are therefore re- 
duced, leaving the period, its significance, the amplitude 
of the variation, its (ffux excess)-color, and finally the light 
curve shape as classification parameters. The light curve 
shape is of particular interest as it can be described math- 
ematically and parameterized in terms of the parameters 
resulting from low order Fourier fits to the data. For clas- 
sical Cepheids the Fourier parameters show a progression 
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Fig. 14. Theoretical expected number densities of LPVs 
according to Renzini (1998) calculated with the small 
bulge decomposition of (Kent) l)l989f) . The field is cen- 
tered on the nucleus of M31, north is on top, east is on 
the left of the diagram. The WeCAPP field for which we 
calculated the number density in Fig.^jis shown by blue 
lines. Contour levels are 1, 2, 3, 4 and 5 xlO^^, 0.1, 0.5, 
1, 5 and 10 sources [arcsec~^]. The red line marks the 
5 X 10~^ contour level. 



with the period of the var iation, echoing the well-known 
'Hertzsprung progression' ijHertzspruel Il924|) of the light 
curve shape with period. For type II Cepheids likewise 
correlations have already been found; the analysis of the 
RV Tauri stars in our sample reveals a correlation between 
different phase parameters. Fourier decomposition of the 
light curves therefore is a powerful method to support 
the classification of Cepheid-like variables and discrimi- 
nate them from other type of variables. On first ground 
our classification scheme is based on the position of the 
stars in the i?-band period-amplitude plane. For certain 
groups of stars we refine and check our classification by 
using the Fourier parameters of the fits to the light curves. 

4.3.1. Period-amplitude-relations 

The amplitudes AFr and AFj of the variation in the R- 
band and /-band are determined as half the difference 
of the maximum to the minimum of the light curves. As 
we calculate the amplitudes after having eliminated data 
points with the largest error bars and biggest flux differ- 
ences this amplitudes act as a lower limit only, but can 
be regarded as a robust measurement of the variation 
amplitude. The amplitudes were transformed in magni- 
tudes using the R- a nd /-band fluxes of Vega taken from 
iBinnev fc Merrifieldl ((1988,'! . Note that these magnitudes 



' I ' I I I I ' I I I I I I I ' ' ' ' I ' ' ' ' J I I 
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Fig. 15. Period-amplitude relation in the /?-band for 
the catalogue sources (eclipsing binaries, irregular, long- 
variation, and miscellaneous variables excluded). The am- 
plitudes were transformed in magnitudes using the /?-band 
fluxes of Vega. Note that these magnitudes reflect the flux 
difference only and are not the real variation magnitude 
of a single star. In this diagram three populations are vis- 
ible. The black lines show the defining relations (from left 
to right) for group I (population I Cepheids: black dots), 
group II (type II Cepheids: blue dots, RV Tauri stars: red 
dots, RV Tauri candidates: open magenta circles, SR vari- 
ables: green dots), and group III (LPVs). For the RV Tauri 
stars and candidates we use the single or fundamental pe- 
riod (minimum to minimum). The sources in group III are 
color coded according to the significance of the /-band pe- 
riod as returned by the Lomb algorithm (low significance: 
black — blue — > red: high significance). 

reflect the flux difference on the frame only, and are not 
the real variation magnitudes of a single star. 

Figure [T^ shows the relation between the logarithm of 
the period and the logarithm of the amplitude AFn (i.e. 
the variation magnitude as measured in the frame) for the 
sources of the final catalogue. Three different groups can 
be distinguished in this diagram: 

— Group I has periods between 1.7 and 21.5 days, 
there seems to be a correlation of the period and 
the amplitude of the variation in a sense that 
bigger periods show larger amplitudes. We require 
-2.5 log(AFi^/Fvcga) < 23.1 - 2.76 (logP^ - 0.55) for 
sources belonging to this group. This relation as well 
as the dividing relation for group II sources is shown 
as black line in Fig. El These stars are most likely 
connected to the disk, as their light curves and the pe- 
riod regime are connected to Cepheids of type I which 
belong to a young population. 

— Group II has periods in the range between 12 and 
about 140 days. Also in this group the period and 
the variation amplitude are correlated, whereas the se- 
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quence for the group II stars lies at fainter magnitudes 
than the one of the group I stars. The defining relations 
for this group are given by —2.5 log(Ai^fl/i^vcga) > 
23.1-2.76 (log Pfl- 0.55) and -2.5 log(A^^K/Fvega) < 
23.0 — 3.5 (logPft — 1.35). Group II stars most likely 
belong to the old spheroid population; RV Tauri stars, 
type II Cepheids and the low period tail of the semi- 
regular (SR) stars are found in this group. 
— The bulk of the detected variables finally belongs to 
group III, variables with periods longer than about 
50 days. For this group no clear correlation between 
the amplitude and the period can be seen. All sources 
with -2.5 log(Afi?/i^vcga) > 23.0-3.5 (log P^- 1.35) 
belong to this group. 

4.3.2. Fourier fits 

We fitted truncated Fourier series of the form 



C + Ao E^L;^ Ai co8{iLu{t - to) + $i 



(13) 



to the i?-band light curves of groups I and II. C defines 
the baseline of the fit, Aq reflects the overall amplitude of 
the light curve, and Ai and <i>i define the amplitudes and 
phases of the different harmonics. The periods entering 
the fit are the periods returned by the Lomb algorithm. A 
fitt ing order of N = 5 l eads to acceptable fits to the data. 

were the first to calculate ampli- 
tude ratios of the form Rij — Ai/Aj and phase differ- 
ences — — i^j of the parameters of the different 
harmonics. Since for classical Cepheids both definitions 
show a progression with the period as well as a correla- 
tion among one another, Fourier analysis proved to be an 
excellent diagnostic to examine the pulsation properties 
of these stars. Fourier analysis therefor e was wi dely used 
in the past to examine Cep heids fe.g.. ISimon fc Moffctt 
l)l986|) . lAlcock et all l)l999|) '). particularly to distinguish 
between fundamental and fir st overton e pulsators (e.g., 
[Antonello fc Poretti (,1986) . Beaiilieu et al.. (1995l) l. 
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Fig. 16. Amplitude ratio i?2i = A2/A1 (upper panel) 
and phase difference $21 = $2 — 2$i (bottom panel) 
determined from the i?-band data for the population I 
Cepheids (group I) plotted against the period. The ratio 
and phase difference drop as the periods approach the res- 
onance P2/P0 ~ 0.5 at about 10 days and rise again after- 
wards. The i?2i value classifies 29 of the detected Cepheids 
as fundamental mode pulsators (closed circles), whereas 2 
are overtone pulsators (crosses). One of the Cepheids falls 
of the fundamental mode sequence i n the phase difference 
diagram. Following iBeaulieu et al ] 1I1995I) we classify this 
source as intermediate type Cepheid (open circle). The 
beat Cepheid candidate for which we used the dominant 
period (presumably the first overtone) in the analysis is 
marked with an asterisk. 



5. Classes of variables 

5.1. Group I and II - Cepheid-like variables 

Groups I and II are populated by Cepheid-like variables 
(population I Cepheids in group I, type II Cepheids and 
RV Tauri stars in group II), and the small period tail of 
semi-regular variable stars, which start to contribute at 
periods from about 16 days onwards. 

5.1.1. Population I Cepheids 

Classical Cepheids are relatively young, intermediate mass 
population I stars. They pulsate in the fundamental mode 
which discriminates them from s-Cepheids, which are be- 
lieved to pulsate in the first overtone. Both groups of 
population I stars have distinctive light curves, the classi- 
cal showing skewed, the s-type showing smooth sinusoidal 
variations. 



We detect 33 population I Cepheids in our sample. To 
check whether our classification is correct and if there are 
first overtone pulsators (s-Cepheids) amongst the detected 
Cepheids we fitted truncated Fourier series (see Sec. 14.3.^ 
to the light curves of group I. Figure [TBI shows the char- 
acteristic progression of the Fourier parameters with the 
period of the light curve for the classical Cepheids: as 
the amplitude ratio i?2i between the first and second har- 
monic drops the corresponding phase difference shows a 
mild rise. The amplitude ratio -R21 declines until a pe- 
riod of about 10 days is reached, and then starts to rise 
again. Due to the suppression of the second harmonic the 
Cepheid light curves in the vicinity of this period look 
quite sinusoidal. Figure (TCI shows, that the the minimum 
of i?2i is connected to a dramatic change of $21 . Generally 
it is believed that this change is connected to the resonance 
between the fundamental mode and the s econd overtone 
P2/P0 « 0.5 fe.g.. . Simon Moffett. (Il986l) 1. 
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Fig. 17. Power spectrum of the beat Cepheid candidate 
in the i?-band. The two peaks in the spectrum correspond 
to periods of 1.73 days and 2.35 days. 



First overtone pulsators can be detected using the 
amphtude ratio R21, as s-Cepheids showing smaller val- 
ues than their f undam ental mode counterparts. Following 
iBeauheu et"al] l)l995|) we demand R21 < 0.3 for P < 3 
days and i?2i <0.2for3<P<5.5 days for the Cepheids 
to be classified as first overtone pulsators. Using this cri- 
terion we identify 2 first overtone pulsators in our sample. 
One source with a R21 value close to the border to the 
s-Cepheids falls of the fundamental mode sequence in the 
phase difference diagram. FoUowing lBeaulieu et al.l 1 1995^ 
we classify this source as intermediate Cepheid. 

One source is clearly separated from the Cepheid re- 
lation in the period-amplitude plane (see Fig. I15() . This 
variable shows Cepheid like variations which are mod- 
ulated by a period of 208 days in R. The modulation 
is most likely due to another variable source inside the 
PSF. Fourier decomposition of the _R-band data shows 
that the source has two excited periods, 1.73 days and 
2.35 days, besides the subdominant modulation of the 
light curve. In Fig. ll7l we show the power spectrum of this 
source. The two periods would classify the source as a beat 
Cepheid, the first one detected in M31. Beat Cepheids are 
a rare sub-class of Cepheids in which two pulsation modes 
are excited simultaneously. The ratio of the two periods 
Pshort/Piong ~ 0.73 makcs it likely that we see a fun- 
damental mode / first overtone (F/IH) pulsator, if the 
Cepheid hypothesis for this source is confirmed. For F/IH 
beat Cepheids the first overtone should be the dominant 
mode, which is fulfilled for our candidate. Interestingly, 
the position in the (Pi/Pq) - log (Pq) diagram would place 
this beat Ceph eid on a seque nce defined by the SMC beat 
Cepheids of Be aulieu et al. 

I §997), well above the rela- 
tions for the LMC and the Galaxy. This would point to 
a Cepheid in M31 of approximately the same metal con- 
tent by mass as the SMC beat Cepheids. The object is 
resolved in our referenc e image and is also classified as 
a star (261262) in the iHaiman et al.l l|l994l) catalogue. 
Additionally it wa s detected in X-rays and appears in the 
iKong et all l)2002|) catalogue (J004301.8-f411726, see also 
Table ID). 

In addition to the beat Cepheid there is another pe- 
culiar source in group I which remains unclassified at this 
stage. It shows two periods of 6.95 days and 20.5 days. 
We will further investigate that object, together with the 



beat Cepheid candidate, and also present a more detailed 
description of the Fourier parameters of group I stars in a 
future publication. 



5.1.2. Type II Cepheids and RV Tauri stars 

Type II Cepheids are low-mass population II variables 
which follow a period luminosity relation about 1.5 mag 
below the classical Cepheid relation. They are found in old 
populations like globular clusters, the halo or the bulge. 
Type II Cepheids with periods less than 5 days are also 
called BL Her or CVB stars, whereas type II Cepheids 
with periods between 10 and 20 days are often referred 
as W Virgini stars. At the upper period limit often a pe- 
riod doubling can be observed in the light curves. An ap- 
proxi mate period limit of about 20 days I Alcock et all 
Il998|) separates type II Cepheids from RV Tauri stars 
which share the same light curve, and o ften also the same 
chemical and dynamical characteristics l)FokinLl200lh . RV 
Tauri stars can be recognized by a typical double-wave 
light curve with alternating deep and shallow minima, 
their semi-periods (minimum to minimum) ranging from 
20 to over 50 days. The current understanding places RV 
Tauri stars at the end of stellar ev olution. After lea.ving th e 
Asymptotic Giant Branch fAGB. Ilben fc Renzinil l(l983l) 'l 
they move left in the HR-diagram entering the instability 
strip at high luminosities. For a recen t review on popula - 
tion II Cepheids and related stars see IWallersteinI l)2002j) . 
A comprehensive collection of type II Cepheids and RV 
Tauri stars detected in the LMC by the MA CHO collab- 
oration can be found in lAlcock et al.l l(l998t) . 

We detect 37 RV Tauri stars and 11 RV Tauri can- 
didates which makes this catalogue to one of the biggest 
collections of RV Tauris to date. The light curves show 
the typical alternation of deep and shallow minima, the 
second maximum being fainter than the first one, al- 
though for a few sources this latter difference is rather 
marginal. The light curve shape can be divided into 
two groups, one resembling the 'flat-topped' Cepheid II 
shape (see Fig. I26II . the other one showing sinusoidal 
variations. To avoid a mis-classification of the sinusoidal 
light curves of RV Tauri stars with the not too differ- 
ent light curves of /3-Lyrae eclipsing binaries, we extract 
the phase difference between the second and the fourth 
h armonic from Fourier fits to the light c urves. According 
to lSzvmanski. Kubiak. fc Udalskil l(200lh these two phases 
should be strongly coupled for sinusoidal contact systems, 
yielding A$ = <I>4 — 2$2 — (M. Szymanski, priv.comm.). 
For pulsating stars the correlation is much weaker. 5 of our 
RV Tauri candidates with sinusoidal light curves have a 
A$ which is within 1-a compatible with (or equivalently 
2tt), another 6 lie within 3-a. We classify all 11 RV Tauris 
whose A$ values are within S-cr compatible with as RV 
Tauri candidates, the remaining 37 sources as RV Tauri 
stars. 

Using the fundamental period (minimum to minimum) 
for the Fourier extension of the RV Tauri sources yields rel- 
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Fig. 18. Phase differences $ij = $i — i (f>j as determined 
from the i?-band data plotted against each other. Black 
dots: population I Cepheids (group f), blue dots: type II 
Cepheids (group II), red dots: RV Tauri stars (group II), 
open magenta circles: RV Tauri candidates (group II). As 
green dots we sh ow the LMC type I I Cepheids and RV 
Tauri stars from lAlcock et al. which fall on the 

sequences of the WeCAPP sources. For the RV Tauri stars 
we use the formal period in the analysis. The classical 
Cepheids show a clear correlation of the two parameters, 
but also for the type II Cepheids and the RV Tauri stars 
both phase differences are correlated. The sequences of RV 
Tauri stars and type II Cepheids overlap, making the RV 
Tauri sequence to an extension of the Cepheid II sequence. 
This favors the close connection between these two types 
of stars. 



atively bad fits and subsequently more uncertain Fourier 
parameters. We therefore use the formal period in the 
analysis which yields a more appropriate description of the 
data, better fits, and more reliable Fourier parameters. In 
Fig- El we show a 'phase-phase' diagram as a further re- 
sult of the Fourier analysis of the RV Tauri and candidate 
RV Tauri light curves. The phase differences $31 and $21 
correlate for the RV Tauri light curves (red dots) show- 
ing that these sources form a homogeneous group of stars. 
Also the RV Tauri candidates (magenta open circles) fol- 
low the relation supporting the RV Tauri nature of these 
stars. Furthermore we show the detected type II Cepheids 
as blue dots in this figure. The sequences of RV Tauri stars 
and type II Cepheids overlap, making the RV Tauri se- 
quence to an extension of the Cepheid sequence at higher 
periods. This supports the close connection between RV 
Tauri stars and type II Cepheids. In fact, lAlcock et alJ 
showed that a single period-luminosity-color rela- 
tionship describes both the type II Cepheids and RV Tauri 



Period [days] 

Fig. 19. Histogram of the periods of the Cepheid vari- 
ables in group I. 



stars in the LMC. We show the phase parameters for the 
stars presented in this study as green dots in Fig. 17. For 
this purpose we re-analyzed the light curves taken from 
the MACHO project database'^ and used the formal period 
to derive the Fourier parameters. Two stars with an un- 
certain classification were rejected from the analysis. The 
LMC variables lie on the sequence of the WeCAPP RV 
Tauri and Cepheid II stars, supporting our classification of 
the variables of group II. We will present a more thorough 
discussion of the Fourier parameters of the Cepheid-like 
variables in group I and II in a future publication. 

The period range of the 93 type II Cepheids extends 
to periods larger th an the approximate limit of about 20 
days as proposed bv lAlcock et alJ 1 19981) . There is still the 
possibility that some of these long period Cepheids are in 
reality RV Tauri stars, since noise in the light curves can 
prevent us to detect the alternation of deep and shallow 
maxima as required for the assignment of the variable as 
a RV Tauri candidate or star. The majority of the light 
curves of this type of variables are of ^flat-topped' shape 
ijKweei Il967l) with relatively long and flat maxima. As 
already mentioned, the light curves of type II Cepheids 
show a progression of the phase differences and occupy 
distinct places in Fourier space. The sample of type II 
Cepheids outnumbers the sample of population I Cepheids 
by about a factor of three. As the type II Cepheids trace 
the old bulge population this is not unexpected, even when 
taking into account their smaller brightness (at the same 
period) if compared to population I Cepheids. The pop- 
ulation I Cepheids, more massive and younger than the 
type II Cepheids, are typical members of the M31 disk 
population, which is superimposed to the M31 bulge. 

The third constituent of the group II sources are small 
period semi-regular variables which contribute from peri- 
ods about 16 days onwards. The Fourier analysis of the 
light curves shows no correlations of the amplitude ratios 
or phase differences with period. Also phase-phase dia- 
grams do not reveal a correlation of the phase parameters. 

In Fig. 1191 we give the period distribution of the popu- 
lation I Cepheids, whereas Fig.l2UI shows the distribution 
of periods for variables belonging to group II. In red, ma- 
genta, blue, and green we present the distribution for the 
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Fig. 20. Histogram of the periods of variables belonging 
to group II. We show the distributions for type II Cepheids 
(blue line), RV Tauri stars (red line), RV Tauri candidates 
(magenta line) and the semi-regular stars (green line). For 
the RV Tauri stars and candidates we use the single or 
fundamental period (minimum to minimum). 

RV Tauri stars, RV Tauri candidates, type II Cepheids, 
and semi-regular variables. 

5.2. Group III - LPVs 

Group III consists of Long Period Variables (LPVs), i.e. 
Mira and semi-regular variable stars. LPVs are members 
of the AGB which marks the final stage of the stellar evo- 
lution for intermediate mass stars with masses between 
0.5 and 8 Mq. These stars evolve from the main sequence 
and populate the red giant branch up to a maximum lu- 
minosity at the TRGB (tip of the red giant branch). After 
ignition of helium core burning they drop in luminosity 
and form the horizontal branch. At the end of this sta- 
dium the luminosity again is rising as the stars evolve up- 
wards the AGB. This evolution sequence is characterized 
by pulsation and extensive mass loss. 

LPVs are a very promising tool for various astrophys- 
ical questions. They follow tight period-lum inosity rela- 
tions in the near-IR (especially the K -band) ijPeast et al.l 
Il989t IWoodl. 120001 IPeast et all l2002|) which makes them 
excellent galactic and extragalactic distance estimators. 
They are furthermore good indicators of the parent pop- 
ulation to which they belong, as there exists a depen- 
dence of the luminosity (hence period) of Miras on the age. 
Longer period Miras should have higher mass progenitors 
and therefore belong to a younger population. LPVs are 
historically separated in two main groups: i) Miras with 
regular variations, periods between 80 and 1000 days and 
an amplitude of the variation in the V^-band of more than 
2.5 mag; ii) semi-regulars (SR) with less regular variations, 
smaller periods and a 1^-band variation smaller than 2.5 
mag. Semi-regulars are dived in two groups, SRas with 
more regular variations and SRbs with variations less reg- 
ular. 

Recently iLebzelter et al.l ll2002l) proposed a new clas- 
sification scheme less dependent on this kind of artificial 
division and used it for the AGAPEROS survey of vari- 
able red stars towards the Magellanic Clouds. Their clas- 
sification is based on the regularity of the variation alone, 
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Fig. 21. Distribution of the periods of LPVs in the bulge 
of M31 detected by the WeCAPP project. Top: class 

i - regular and semi-regular variables (significance level 
P(> z) < 10"^" in i? or /). Bottom: class ii - irregular 
variables (significance level P(> z) > 10"^" ini? and /). 
Only irregulars selected by the significance cut are shown. 
Note that the distribution for the sources classified as ir- 
regular variables is of qualitative nature only as the peri- 
ods for these sources are not well determined. In general 
irregular variables show smaller 'periods' than the group 
of regular and semi-regular variables. 

providing three classes of stars: i) LPVs with regular vari- 
ation ii) LPVs with semi-regular variation iii) LPVs with 
irregular variation. No cuts in amplitude were applied, 
therefore class i contains members of the classical Mira 
group as well as members of the SRa type of objects. 

We modify this classification scheme for our purposes 
and define a significance cut of 10^^° according to the 
period finding algorithm in i? or / for the division of our 
sample of LPVs in classes i (regular and semi-regular) and 

ii (irregular). Class ii coincides with the class of irregular 
stars introduced in Sec. 0| Because of the non-linear de- 
pendence of the significance on the S/N of the light curve 
(see also discussion below) we abandon to invent an au- 
tomatic cut for the division into semi-regular and regular 
variations. The cut at 10"^" is somewhat subjective as also 
stars with semi-regular, but low S/N light curves are clas- 
sified as irregular, as the low S/N prevents the period to 
be determined with better significance. On the other hand 
this cut ensures that most of the irregular light curves are 
classified in the right way. The histograms of the periods 
for the irregular and the regular/semi-regular sample of 
LPVs are shown in Fig. [21 The irregular variables show 
in the mean smaller 'periods' as the group of regular or 
semi-regular variables. Note, that this statement is of qual- 
itative nature only, as the periods for these sources are not 
well determined. 

We have color coded the LPVs in Fig. El according 
to the significance of the /-band period. The slope of the 
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at 




R-l [m'agfj 

Fig. 22. Variation color R — I shown as a function of the 
i?-band variation magnitude. The sources of groups I and 
II are shown as big closed circles, the sources of group III 
as small closed circles. The color coding is the same as in 
Fig. El i-e. the LPVs are again color coded according to 
the significance of the /-band period. The sources in group 
I and II as a whole show bluer variations than the LPVs 
in group III. The significance of the /-band period rises 
for brighter R variation magnitudes and redder variation 
colors. This is a result of the higher S/N of the /-band 
light curves, but can in part also be attributed to the 
enhanced regularity of the light curves. 




log(Period) [days] 

Fig. 23. This figure shows the variation color /? — / as 
function of the period of the variable sources. The color 
coding is the same as in Figs. 1151 and 1221 Note, that due 
to the normalization to the fluxes of Vega in the par- 
ticular system this R — I color has not to be mistaken 
with the color of the variation amplitudes of the sources. 
Nevertheless, some trends are visible in this diagram: As 
the variation of the LPVs gets redder with increasing pe- 
riod, the /-band periods become more significant, reflect- 
ing the rising S/N in I due to higher R — I values and a 
rising regularity of the light curves. 



variation magnitude-period relation changes from slightly 
positive for the low significance tail (black dots) to nega- 
tive for the sources with intermediate and high significant 
periods (green and red dots). The same trend can be seen 
in Fig. |52 which shows the variation color i? — / as a func- 
tion of the /?-band variation magnitude. The significance 
of the /-band period rises for brighter R variation magni- 
tudes and redder variation colors. The differences in the 
significance of the LPV /-band periods therefore is a re- 
sult of the different S/N of the /-band light curves (rising 
S/N in / due to higher R — I values), but can in part be 
also attributed to different stages of regularity in the light 
curves. Figure finally shows the variation color R — I 
as function of the period of the variable sources. The vari- 
ation of the LPVs gets redder with increasing period, and 
at the same time shows more significant periods. 

6. The Catalogue 

The final catalogue comprises 23781 entries with group I 
containing one beat Cepheid candidate, 2 s-Cepheids, one 
Cepheid of intermediate type, and 29 classical Cepheids 
pulsating in the fundamental mode. One further group I 
source which remains unclassified in this catalogue shows 
two periods. 

Group II comprises 93 type II Cepheids, 37 RV Tauri 
stars, 11 RV Tauri candidates, and 193 low period semi- 



regular variables. Group III consists of 4287 irregular 
and 18974 regular /semi-regular variables. 82 presumably 
group III members show variations on timescales longer 
than the survey length. The positions in the i?-band 
period- amplitude plane of the 31 eclipsing binary candi- 
dates coincide with the positions of groups I and II. Finally 
we detect 39 miscellaneous variables, among them 16 no- 
vae and 15 R Coronae Borealis candidates. 

The full catalog is available in electronic form at CDS. 
We give the name, the position in the WCS, the periods 
(if available) derived from the R- and /-band data, the 
amplitudes of the variation AFn and A//, and provide 
a classification according to Sees. 0] and El DC (classical 
or 5-Cepheids), SC (s-Cepheids) , BC (beat Cepheids), IC 
(intermediate Cepheids), W (type II Cepheids), RV (RV 
Tauri stars), S (regular or and semi-regular variations), 
and I (irregular variables). We also mark the identified 
Novae (N), eclipsing binary candidates (E), RGB candi- 
dates (RGB) and other miscellaneous variables (M). As an 
illustration of its content we list in Table|21the entries 100 
— 120 in the catalogue. In Figs. [^to|221we show typical 
light curves for each of the groups. 

7. Correlation with other catalogues 

We cross-correlate our /?— band selected catalogue 
with the General Catalogue of Variable Stars (GCVS, 
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Id 


q(2000.0) 


5(2000.0) 


Pr [days] 


-log(sfl) 


Pr [days] 


-log(S7) 


AFr 


AFi 


Type 


WeCAPP.VOOlOO 


0h42ml8.46s 


41d08'09.0" 


119.70 


13.18 


122.10 


8.52 


0.37 


1.06 


I 


WeCAPP.VOOlOl 


0h42m09.65s 


41d08'08.1" 


274.99 


21.65 


277.50 


26.10 


0.78 


2.86 


S 


WeCAPP_V00102 


0h42m3 1.76s 


41d08'11.5" 


500.40 


9.60 


478.70 


22.69 


0.47 


3.48 


s 


WeCAPP_V00103 


0h42m20.06s 


41d08'04.4" 


180.79 


29.32 


180.21 


26.71 


0.76 


1.70 


s 


WeCAPP_V00104 


0h42m31.53s 


41d08'08.l" 


289.96 


18.52 


293.51 


23.71 


0.74 


3.18 


s 


WeCAPP_V00105 


0h42m52.12s 


41d08'11.0" 


299.24 


33.60 


302.18 


41.93 


0.63 


2.59 


s 


WeCAPP_V00106 


0h43ml7.73s 


41d08'l2.8" 


136.3 


8.920 


140.74 


4.86 


0.49 


1.05 


I 


WeCAPP_V00107 


0h42ml9.93s 


41d08'09.7" 


284.59 


10.99 


273.54 


29.43 


0.43 


2.03 


s 


WeCAPP_V00108 


0h42m45.51s 


41d08'10.7" 


268.01 


19.64 


267.94 


38.87 


0.57 


2.56 


s 


WeCAPP_V00109 


0h42ml8.84s 


41d08'08.8" 


240.35 


21.15 


247.07 


27.50 


0.51 


1.74 


s 


WeCAPP.VOOllO 


0h42ml7.75s 


41d08'08.8" 


423.96 


5.730 


423.96 


24.32 


0.48 


1.36 


s 


WeCAPP.VOOlll 


0M2ml7.43s 


41d08'08.6" 


280.04 


16.27 


278.51 


27.95 


0.48 


1.69 


s 


WeCAPP_V00112 


0h42m30.69s 


41d08'07.9" 


216.48 


17.07 


218.03 


21.92 


0.81 


2.21 


s 


WeCAPP_V00113 


0h42ml7.04s 


41d08'06.8" 


229.50 


25.36 


229.51 


28.31 


0.79 


2.58 


s 


WeCAPP_V00114 


0M3m02.66s 


41d08'13.0" 


226.44 


29.33 


230.99 


35.27 


0.53 


1.96 


s 


WeCAPP_V00115 


0h42m35.32s 


41d08'12.0" 














L 


WeCAPP_V00116 


0h42m23.24s 


41d08'll.l" 


314.68 


24.78 


319.13 


31.99 


0.53 


2.60 


s 


WeCAPP_V00117 


0h42ml4.70s 


41d08'l0.l" 


214.96 


26.14 


209.84 


27.26 


0.59 


1.77 


s 


WeCAPP_V00118 


0h42m24.75s 


41d08'12.4" 


277.49 


10.20 


280.04 


23.22 


0.42 


1.53 


s 


WeCAPP_V00119 


0h42m47.43s 


41d08'll.0" 


497.12 


24.86 


505.28 


28.47 


0.80 


6.19 


s 


WeCAPP_V00120 


0h42m24.29s 


41d08'll.5" 


198.21 


21.46 


198.61 


23.74 


0.46 


2.82 


s 



Table 2. Extract from the WeCAPP catalogue of variable stars. We give the identification, the right ascension and 
declination, the assigned period and its significance s (0.1 = un-significant, 70 = highly significant) in R and /, the 
amplitude of variation in both bands, and finally a classification according to the scheme presented in Sees. 0] and |5| 
The variation amplitudes AFr and AFj are given in units of [10^^ Jy]. 
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Fig. 24. Light curve of a (5-Cepheid (group 1) in the R- 
band (top panels) and the /-band (bottom panels). In the 
upper panels we show the light curve i?-band, the lower 
panels show the light curves in the /-band. In both bands 
the light curves are presented in the time domain (top 
for each band) and and in the phase domain, i.e. con- 
volved with the derived periods (bottom for each band). 
Additionally we show the periods Pr and Pj derived in 
the R- and /-band, respectively. 

IPurlevich et al.' (l99<f}) and catalogues selected in the X- 
ray (Chandra: .Kaaret. (,20021) fc .Kong et al.. (20021 ). 



P„=29 06 days 





loce iiiK) laiQi 




Fig. 25. Light curve of a long period type II Cepheid 
(group II) in the /?-band (top panels) and the /-band (bot- 
tom panels). 

We find 23 coincidences with the iKaarell l)2002|) cata- 
logue when using a search radius of 1". To estimate how 
many coincidences we expect by chance we follo w the ap- 
proach proposed bv lHornschemeier et al.l l)200l|) . We shift 
one of the catalogues to be compared (in our case we shift 
the catalogues taken from the literature) by 20" in north- 
east, south-east, south-west and north-west direction, and 
check for coincidences between these mock catalogues and 
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WeCAPP 


q(2000) 


5(2000) 


Kaaret 


2MASS 


A r(arcsec) 


10578 


0h42ml0.30s 


41dl5'l0.4" 


J004210.2+411510 


0042102+411510 


0.5 


16262 


0h42ml2.14s 


41dl7'58.8" 


J004212. 1+411758 


0042121+411758 


0.4 


8624 


0h42ml8.68s 


41dl4'01.8" 


J004218.6+411402 


0042186+411402 


0.5 


8946 


0h42m21.58s 


41dl4'19.8" 


J004221. 5+411419 


0042215+411419 


0.4 


7582 


0h42m25.11s 


41dl3'40.7" 


J004225. 1+411340 




0.6 


19068 


0h42m31.26s 


41dl9'38.7" 


J004231. 2+411938 


0042312+411938 


0.2 


7979 


0h42m36.59s 


41dl3'49.9" 


J004236.6+411350 




0.7 


9082 


0h42m39.60s 


41dl4'28.6" 


J004239.5+411428 




0.9 


11406 


0h42m41.45s 


41dl5'23.9" 


J004241. 4+411524 




0.3 


16322 


0h42m43.97s 


41dl7'55.5" 


J004243.9+411755 




0.4 


15344 


0h42m44.88s 


41dl7'39.5" 


J004244.8+411740 




0.4 


12503 


0h42m46.99s 


41dl6'15.3" 


J004246.9+411615 




0.7 


8501 


0h42m47.13s 


41dl4'13.9" 


J004247.1+411413 




0.8 


10431 


0h42m47.45s 


41dl5'07.6" 


J004247.4+411507 




0.2 


11736 


0h42m47.90s 


41dl5'50.6" 


J004247.8+411550 




0.6 


12891 


0h42m48.65s 


41dl6'25.0" 


J004248.6+411624 




0.7 


20170 


0h42m55.27s 


41d20'45.l" 


J004255. 3+412045 




0.7 


18374 


0h42m59.66s 


41dl9'19.3" 


J004259.6+411919 




0.3 


12404 


0h42m59.88s 


41dl6'06.0" 


J004259.8+411606 




0.1 


21059 


0h43m03.31s 


41d2l'21.8" 


J004303.2+412121 




0.5 


6049 


0h43m08.42s 


41dl2'46.9" 


J004308.4+411247 




0.7 


8360 


0h43m09.88s 


41dl9'00.8" 


J004309.8+411900 




0.7 


9991 


0h43ml0.62s 


41dl4'51.3" 


J004310.5+411451 


0043106+411451 


0.5 



Table 3. Coincidences between the WeCAPP catalogue of variable stars and the X-ray selected catalogue of point 
sources bv lKaare3 1|2002() . A sear ch radius of 1" was used. We give the identifier and WCS coordinates for t he WeCAPP 
sources, and the identifier for the' Kaare^ l)2002(l sources. For sources coincident with the 2MASS catalogue ijCutri et all 
we give the 2MASS identifier also. Additionally we show for each of the correlated sources the difference A r 
of the matching of both catalogues. S ource WeCAPP V104 31 is identical with the (dwarf-)nova WeCAPP-N2000-05 
which shows supersoft X-ray emission ((Pietsch et al.L 120051) . 




Fig. 26. Light curve of a RV Tauri star (group II) in the 
i?-band (top panels) and the /-band (bottom panels). We 
show the light curves convolved with the formal periods 
(deep minimum to deep minimum) . The double wave light 
curve with alternating deep and shallow minima is nicely 
uncovered. 



Fig. 27. Light curve of a low period semi-regular vari- 
able (group II) in the i?-band (top panels) and the /-band 
(bottom panels). 



our catalogue. By averaging these numbers we get the ex- 
pected number of coincidences by chance. 

Applying a search radius of 1" the number of false co- 
incidences with the iKaarei] 1 20021) catalogue becomes 12. 
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Pp,= 283.19 days 




1- e 



P,=2S5,S2 days 




Fig. 28. Light curve of a LPV with regular variation 
(group III) in the i?-band (top panels) and the /-band 
(bottom panels). As illustration of the agreement of the 
data taken with different telescopes we show the data 
collected at Calar Alto as red dots, and data taken at 
Wendelstein as black dots. 




Fig. 29. Light curve of an irregular LPV (group III). Top 
panel: i?-band. Bottom panel: /-band. 



This high level of false coincidences (about 50%) results 
from the high number density of the variable star cata- 
logue. Nevertheless, the false detections rate suggests that 
about 11 of the coincidences should be real. This is sup- 
ported by the fact, that 6 of the co incident sources h ave 
an entry in the 2 MASS catalogue l)Cutri et alJ.l200(]() . 8 
of the Kaaret- WeCAPP sources are ident ified a s globular 
clusters (4 of them with a 2MASS entrv. iKaaret (2002)). 
Furthermore we find tha t 12 of the s ource s have coinci- 
dent counterparts in the 'Kon g et alJ (jifloij) catalogue as 
well (see below). In Table El we give the positions of all 23 



Fig. 30. Light curve of a LPV (group III) with changing 
variation amplitudes from cycle to cycle. Top panels: /?- 
band. Bottom panels: /-band. 




Fig. 31. Light curve of a R Coronae Borealis candidate 
(miscellaneous variables). Top panel: R-hwA. Bottom 
panel: /-band. These rare carbon-rich stars show unpre- 
dictable and frequent fading of light in optical wavebands, 
most probably due to dust grains ejected from t he stellar 
photo sphere. For a review on these RGB stars see lClavtonl 



coincidences, for sources with entries in the 2MASS cata- 
logue we give the 2MASS name also. Two of the Kaaret- 
WeCAPP sources with 2MASS fluxes (WeCAPP_V8946 
and We CAPP_V10578 ) have no globular cluster (gc) coun- 
terpart l|Kaaretll200^ . Whereas V8946 coincides with a 
gc-candidate identified bv lWirth. Smarr. fc Brunei l)l985f) 
a gc-counterpart for WeCAPP_V10578 is completely un- 
known. 
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WeCAPP 


a(2000) 


5(2000) 


Kon 




A r(arcsec) 


14906 






41 HI 7'9n 


2" 


7004207 


04-41 17iq 


n 

u 


q 

y 


15982 




49s 


41 HI 


u 


1004209 


4+411745 



u 


7 


16262 




14s 


41 HI 7' ^8 


8" 
o 


1004212 


1 -i_4i 1 758 


u 


a 
o 


5801 






41 HI 9''^4 


8" 
o 


7004215 


n+41 19'^4 



u 


u 


14846 






41dl7'21 


q" 


1004215 


6+411721 


n 

u 


7 


12198 






41H1 ^'^"^ 


o 


1004216 


n+41 1 559 


n 

u 


4 


7582 




lis 


41 HI 'H'4n 


7" 


T004225 


1 _|_41 1 S40 


n 

u 


Q 
O 


20282 




61s 


41 H9n'4Q 


n" 

u 


1004227 


6+412048 



u 


Q 
O 


13905 




23s 


41 HI fi''i4 


n/l 
O 


7004230 


9_1_41 



u 


q 

y 


19068 




26s 


41 HI Q''^8 


7" 


7004231 


9_i_4i iq'^q 



u 


7 


11863 




51s 


41H1 ^Mii 


7" 


J004232 


5+411545 


n 

u 


2 


10895 






41 HI 'i'99 


q" 

y 


T004234 


74-41 1 523 


n 

u 


Q 
O 


7979 






41 H 1 V4Q 


q" 
y 


1004236 


54-41 1 350 


n 

u 


Q 
O 


11406 




45s 


41H1 ^'9"^ 


q" 

y 


7004241 


4+41 1523 



u 


Q 
O 


10626 




74s 


41H1 ^'14 


4" 


7004243 


7+411514 



u 


a 
o 


15344 






41 H 1 7''^Q 


u 


7004244 


8+41 17'^q 


n 
u 


■J 
o 


15318 




12s 


41 HI 7''^(S 


nil 
O 


7004246 


n+41 17'^fi 


n 

u 


8 
o 


12503 




QQs 


41 HI fi'l 


o// 
O 


7004246 


9+41 1615 


n 

u 


p. 
u 


11736 






41H1 ^'"Ml 


u 


7004247 


8+411549 



u 


8 
o 


22406 




94s 


41 H94'nfi 


1" 


7004248 


q-l_41 2406 


I 


n 

u 


17309 






41 HI 8'44 


4" 


7004256 


q-l_41 1 844 



u 


2 


18374 






41 HI Q'1 Q 


O 


7004259 


6+411919 


n 

u 


4 


12404 




SSs 


41H1fi'n(S 


n" 

u 


7004259 


8+41 1 606 


n 

u 


2 


15035 




sris 


41 HI 7'9fl 


u 


7004301 


8+41 1726 


n 

u 


5 


10871 
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41 HI 'i'99 


o// 

o 


7004302 


9+411522 



u 


p. 
u 


21059 
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/1 1 H91 '91 
4iQzr zi 


8^' 
o 


7004303 


•^-l_41 21 22 


n 
U 
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0h43m03 
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41dl8'04 


7" 


7004*^0'^ 


On^'y: _L _Louu 





5 


19890 


01i43m07 


55s 


41d20'20 


0" 


J004307 


5+412020 





1 


17815 


01i43m09 


88s 


41dl9'00 


8" 


J004309 


8+411901 





4 


9991 


01i43ml0 


62s 


41dl4'51 


3" 


J004310 


6+411451 





8 


17140 


01i43ml6 


14s 


41dl8'41 


5" 


J004316 


1+411841 





6 



Table 4. Coincidences betwe en the WeCAPP catalogue of variable stars and the X-ray selected catalogue of point 
sources by iKong et al ] lj2nn2|> . A search radius of I" was u s ed. W e give the identifier and WCS coordinates for 
the WeCAPP sources, and the identifier for the iKong et"al sources. Additionally we show for each of the 

correlated sources the difference A r of the matching of both catalogues. The position of the beat Cepheid candidate 
WeCAPP_15035 coincides with the X-ray source J004301.8+411726. 



The iKong et al catalogue contains 31 coinci- 

dences with our variable star catalogue when applying a 
1" search radius. As we expect 17 coincidences by chance 
about 14 of the coincidences should be real. As me ntioned 
above 12 of the sources are also coincident with the lKaareJ 
l)2002() catalogue. We show the position of all 31 coinci- 
dences in Table 0] 

The GCVS l)Durlevich et all Il996(l comprises 250 en- 
tries within our field of view, most of them being of nova 
type. 161 (6) were classified as novae, additionally 53 as 
novae of fast type (NA), 9 (2) of slow type (NB), and (2) 
of very slow type (NC). Furthermore the GCVS reports 1 
SNIa remnant in our field, 1 (1) LC type (irregular vari- 
able supergiant), 3 (1) semi-regular variables (SR), 3 S- 
Cepheids, (1) SDOR variable, 2 (1+2) irregular variables, 
and 1 unstudied (S: classification) variable. The numbers 
in brackets give additional candidates for the particular 
variable star group with questionable classifications ac- 
cording to the GCVS. 



The crowding of the WeCAPP sources demands that 
we keep the search radius at 1", even when comparing 
with the GCVS, whose accuracy is well below this value. 
We find 27 coincidences, the number of false detections 
(21) suggests however, that most of the coincidences are 
not real. As furthermore 23 of them are classified as novae 
in the GCVS we only show the 4 coincidences which we 
regard as real in Table [S] Note, that we do not match 
the positions of all three (5-Cepheids found in the GCVS 
inside our field. The GCVS Cepheid V0934 has a position 
difference of 2.1" when compared to our position, GCVS 
Cepheid V0811 is off by 10.2". Because of these in parts 
large uncertainties in the positions of the GCVS we are 
not trying to match the two remaining semi-regulars given 
their smaller variation amplitude. 

8. Conclusions and outlook 

We presented the catalogue of variable stars detected by 
the WeCAPP project in the central parts of M31. The ob- 
servations in the optical R and / bands towards the bulge 
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WeCAPP 


a(2000) 


5(2000) 


GOVS 


A r(arcsec) 


GOVS class 


WeCAPP class 


906 


0h43m22.50s 


41d07'28.6" 


V0871 


0.7 


I: 


I 


2375 


0h43ml5.73s 


41d07'54.7" 


V0727 


0.1 


lA: 


I 


559 


0M3m24.16s 


41d07'12.7" 


V0743 


0.2 


DCEP 


DC 


742 


0h43ml7.27s 


44dl5'58.0" 


V0726 


1.0 


SR 


I 



Table 5. Co incidences between th e WeCAPP catalogue of variable stars and the General Catalogue of Variable 
Stars rCCVS. Irwievich et a,lJ l|l99(ill '). A search radius of 1" was used. We give the identifier, WCS coordinates and 
a classification for the WeCAPP sources, and the identifier and a classification according to the GCVS. Additionally 
we show for each of the correlated sources the difference A r of the matching of both catalogues. 




the LMC bar, i s populated by ell ipsoidal red giants 
i Soszvnski et alll2004 . IWoodI |l999l) found five distinct 



I9I» 



Fig. 32. Light curve of a nova (miscellaneous variables). 
Top panel: i?-band. Bottom panel: /-band. A paper con- 
cerning the characteristics of all detected novae is in prepa- 
ration. 

of M31 covering the years 2000-2003 with very good time 
sampling resulted in a database of over 23000 variable 
sources. The catalogue containing regular, semi-regular 
and irregular LPVs, Cepheids, RV Tauri stars, eclipsing 
binary and R Coronae Borealis candidates enlightens the 
rich population of variable sources in M31. Two more pa- 
pers dealing with the Fourier parameters for the Cepheid- 
like stars and the WeCAPP nova catalogue are in prepa- 
ration. A fraction of the nova catalogue has already been 
correlated with supersoft X-ray sources (SSSs) in M31, 
showing that classical nov ae constitute the major class of 
SSSs in M31 jPietsch et al.. 20051. 

The distribution of the identified variable sources 
shows an asymmetry due to the enhanced extinction in 
the spiral arms projected on the M31 bulge. Assuming 
that the non- variable stars in the bulge behave in the same 
way as the variable counterparts this asymmetry enlarges 
the expected asymmetry in the microlensing signal result- 
ing from the high inclination of M31. Therefore theoret- 
ical calculations have to take a modified distribution of 
sources (or the distribution of the extinction present in a 
given band) into account. 

Recent ly it w a s pro posed that sequence E, first de- 
tected by IWoodI l|l999l) in the MACHO data towards 



period-luminosity sequences, of which three could be at- 
tributed to different evolutionary phases and pulsation 
modes. For two of them, the sequences D and E, it was 
suggested that they are populated by contact binaries, 
and semi-detached binaries with an invisible companion. 
Analyzing the data of the LMC ellipsoidal candidates from 
the OGLE web archive"* and accounting for the distance 
modulus of M31 proves that these sources are below the 
detection limit of our survey. Therefore, we can not put 
any constraints on a possible population of ellipsoidal red 
giants in M31. 

The more detailed examination of the group 111 LPVs 
will be subject to future work. Until now we measured 
the dominant period in these stars which are known to 
show multi-periodic phenomena. Future work will include 
the extraction of the sub-dominant periods and the distri- 
bution of the ratio of the different periods. Furthermore, 
the distribution and possible correlations of the Fourier 
parameters for these stars is an interesting and until now 
unstudied problem which we can address with this dataset. 
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Appendix A: Survey window function and leakage 
function 

Effects induced by the discrete time sampling of astro- 
nomical surveys can yield strong power in frequencies not 
caused by a periodic or quasi-periodic signal. The frequen- 
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cies affected are often connected to the daily separation 
of observations or to observing blocks interrupted yearly 
by the non-observability of the object. In case of M31 the 
observation gap usually lasts from the mid of March to 
the beginning of June. 

A quantitative way to describe the effects of the time 
sampling of the observations on the power spectrum is 
given by the window function W{v). W{i') is defined as 
the Fourier transform of the sampling function s{t), where 
s{ti) is of constant value in case of observations taking 
place at times ti and zero otherwise: 



Wiiy) = / s(t)e 



-2TTil/t 



dt 



(A.l) 



In Figs. Em and ESI 

we show the window function 
W(y) (in the frequency and period representation, respec- 
tively) calculated with the Lomb algorithm for the sam- 
pling given in Field Fl. For this purpose we created a high 
(order of 10"'^'') signal-to- noise-ratio {S/N) light curve with 
s{ti) = l + a(ti) and a minimal noise contribution cr(tj) to 
allow for the calculation of the variance of the data (see 
Eqs.HlandOJ). 

Clearly visible is the strong power in periods around 
450 days (at ly = 2.2 x 10""^ days~^), whereas at small pe- 
riods (high frequencies) no 'aliasing' problem is expected. 
The full width half maximum t'FWHM of the main peak 
of the window function W^v) determines the theoretical 
error in the period calculation 



AP; 



(A.2) 



with Aiy = vfwum- According to lRober'te 
Av can be approximated by Av w 1/T with T the total 
time span of the observations if the sampling is not to 
non-uniform. The measured j^fwhm = 0.0013 day"-'^ is 
1.60 times larger than this theoretical value which can be 
explained by the non-uniform sampling of the observations 
(see Fig. 12). With this value Eg. lA. 21 becomes 



AP K. 6.5 X 



100 d 



(A.3) 



An alternative way to illustrate the sampling effect is 
given by the 'leakage function' for different periodic sig- 
nals. The function describes how a signal which is located 
at zero offset 'leaks' into neighboring frequencies. Figure 
lA. 31 shows the leakage function for a periodic signal of the 
form 



-[^Qsin(27r(ti-io)/P) 



(A.4) 



calculated with the Lomb algorithm for the sampling s{t^ 
of field Fl and different periods P . The problem of high 
amplitudes in the power spectrum for frequencies corre- 
sponding to periods in the range between 350 and 550 
days is apparent in this figure also. 

In Fief. lA.4| wc show the leakage function in the period 
representation. For periods between about 350 and 550 




Fig. A.l. This figure shows the window function W(y) 
present in the survey, the amplitude is normalized to 1. 
W{v) is defined as the Fourier transform of the sampling 
function s(t), W(y) = J s(t) e-'^''''''*dt. The main peak 
of the window function has a full width half maximum 
i^FWHM = 0.0013 day~^ which is 1.60 times larger than 
what one w'ould expect fro m the total duration of the ob- 
servations ([Roberts et allll9871 



-I 




'iBB 5IH1 

Fig. A.2. Window function W{i') shown in the period 
representation. Clearly visible is the sampling induced 
strong power in periods around 450 days which give rise 
to the aliasing problems (see also Figs. lA.Sl and lX^ . 



days the power in the second peak in the Fourier spec- 
trum can exceed the power in the first peak which means 
that the derived period is offset from the real one by the 
position of the second peak. 

Appendix B: Monte Carlo simulation 

Additionally we tested the accuracy of the period de- 
termination with Monte Carlo simulations of the survey. 
We look for the dependence of the accuracy for differ- 
ent mean magnitudes and amplitudes at a fixed position, 
hence fixed noise contributions of the galaxy and the sky. 
We choose the mean magnitude Mu of the variable to be 
Mjf = —1,0,1,2 mag and shift them to the distance of 
M31. For each of these mean magnitudes we change the 
amplitude Aji of the variation from 0.5 mag to 2 mag. 
The background (surface brightness of the galaxy -I- sky 
contribution) is set to M^g^fj — 19 mag per area of the 
PSF. In each run we simulate 10^ light curves according 
to 



FniU) = 10 



-0.4(AfH + (mH-AfH)+AHSin(27r(t,-to)/Pi„)) 



(B.l) 
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Fig. A. 3. Leakage function for the sampling given in the 
survey for different periodic signals. The function shows 
how a signal which is located at zero offset leaks into neigh- 
boring frequencies. We calculated this function for a signal 
according to Eq. (|A.4|I and periods of 5, 10, 20, 50, 100, 
200, 250, 300 (cyan curve), 400, 450, 500, 550, 600, and 
700 days (black curves, if not mentioned otherwise). The 
power in the second peak is high for periods of 400 (blue 
curve), 450 (magenta curve), and 500 days (red curve). In 
this period regime the power in the second peak can exceed 
the power in the first peak (at zero frequency difference) 
resulting in an aliasing of periods. The aliasing problem 
gets lower for periods of 550 (green curve) and 600 days 
(yellow curve), and vanishes for even larger periods. The 
aliasing problem is also evident in the Monte-Carlo simu- 
lation of the survey. 




Fig. A. 4. Leakage function in the period representation 
for the sampling given in the survey for different periodic 
signals. Color coding as in Fig lA.3l For periods between 
350 and about 550 days the power in the second peak 
can exceed the power in the first peak (at zero period 
difference) resulting in a possible aliasing of periods. The 
resulting displacement of the true period can be read from 
the X-axis. 

The epochs ti are given by the sampling of field Fl 
(see Fig. El, periods Pin are chosen randomly between 3 
and 700 days, and to, the epoch of mean luminosity is 
randomly distributed between the beginning and the end 
of the survey. Note, that this also ensures that the phase 
is randomly chosen. For each simulated light curve the 
period is derived according to Eq.^ Figure ^Q] shows the 
derived periods Pout as a function of the input period. It is 
clearly seen that the accuracy of the period determination 
decreases at around 500 days, independently of the S/N 
of the light curves. This can partly be explained by the 
sparse time sampling in the last two campaigns. On the 



other hand also aliasing effects due to the yearly observing 
window seem to play a role for periods between 350 days 
and 550 days (see Figs. EH El and IE}. This holds 
for S'/A^-ratios above a certain threshold. Noise present 
in the light curves tends to level out the different peaks 
in the power spectrum which results in a generally bad 
performance of the period determination for light curves 
with low S/N. For higher S/N light curves the accuracy of 
the period determination is good apart from the problem 
that periods around 400 and 500 days tend to be shifted 
to much smaller (i.e. around half of the true) values. 
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Fig. B.l. Result of simulations for decreasing S/N of the 
light curves from the upper left to the lower right corner. 
The mean magnitude of the source changes in steps of one 
magnitude from Mn = — 1 mag in the upper left corner 
to AIr = 2 mag in the lower right corner. Red dots show 
the accuracy of the period determination for light curves 
with an amplitude An = 2.0 mag, blue dots correspond to 
All = 1-0 niag, and black dots have An = 0.5 mag. The 
apparent horizontal lines at large periods result from the 
finite frequency resolution of the period finding algorithm. 
Above a threshold in S/N the accuracy of the period de- 
termination seems independent of the particular S/N of 
the light curves. In this regime only periods larger than 
500 days have a large false detection probability. Below 
the threshold the accuracy breaks down very rapidly, in- 
dependent of the period. In magenta we show the theoret- 
ical expected errors of the period determination according 
to Eq. 



